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Interim Report 


This report was prepared by the University of Alabama in 
Huntsville under Contract Number NAS8-30216, Study to Define 
Points of Entry for Potential Contaminants into Limestone 
Aquifers, for the George C. Marshall Space Flight Center of 
the National Aeronautics and Space Administration. 
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Interim Report 

Determination of Vo Into of Entry for Potential 
Contaminants Into Limestone Aquifers Using 

Thermal Infrared Imagery 


1.0 Introduction 

1.1 Purpose and Approach of Investigation 

The research reported here is being conducted in order to 
find applications of remote -sensing techniques to ground water hydrology. 
Because such an attempt involves the subsurface, an indirect approach 
must be used. Hydrologic applications of remotely-sensed data are 
straight forward and direct for surface waters. Flooded areas, sediment 
plumes, and thermal differences in water bodies are amenable to direct 
sensing. Ground water, on the other hand, is usually several feet to 
several 10 's of feet below the surface and not directly recordable. 

Soil moisture, however, can be related to ground water, and, furthermore, 
under certain circumstances it can affect surface temperature differen- 
tials. The effect of soil moisture on surface temperature and the 
relationship between soil moisture and ground water is exploited in this 
research. 

1.2 Personnel and Materials 

Thermal differences at the earth's surface were recorded 
before dawn on 26 and 27 June 1974 over the southwestern part of Madison 
County and adjacent Limestone and Morgan Counties by a Bendix LN 2 thermal 
scanner mounted on a C-45 airborn platform. Contracting Office's Repre- 
sentative is Mr. Herman Hamby. Mr. Rex Morton was chief of operations; 

Mr. William Scarbrough operated the instruments. Interpretation of the 
imagery was made by Dr. F. L. Doyle. Research assistants were 0. E. 
McCartney and Steven C. Vavra. 

2.0 Summary of Previous Research (Phase One) 

This research continues that started in June 1973 (phase one) and 
reported in November 1973. During that phase of the research natural 
linear features were identified on ERTS imagery of north Alabama, and 



on high altitude daytime thermal Imagery. Tilly it/Uigery way alyO enhanced 
by density slicing on a Data Colot Model 703-32 In order to minimize 
subjectivity in interpretation. 

The second phase of research is based on the lincations recognized 
in phase one. The results of the second phase are described in this 
report . 

3.0 Relations of Present Investigation to other studies. 

3.1 Lineations (lineaments) have been identified by Drahovzal 
and others (1974) in N/.SA-sponsored research involving the application 
of ERTS imagery to geologic and hydrologic problems. Their interpreta- 
tion of the southwest Madison County area is substantially in accord 
with the lineations reported in the first interim report of F. L. Doyle 
(dated 21 November 1973). 

A major linear trend has been named the Anniston lineament 
by Drahovzal and his colleagues. They also report that in general the 
bearings of surface fractures (joints) do not parallel the bearings of 
the lineations (p. 93). This finding re-enforces the belief expressed 
in the first interim report that the lineations are probably associated 
with deep-seated, "basement" faults below the sedimentary veneer. 
Drahovzal and others also report (p. 452) that gravity anomalies they 
have determined in western Madison County can be interpreted as being 
caused by faulting at the stratigraphic level of the Knox Dolomite 
(roughly 1500 feet below the surface). Such a finding is also in 
accord with the suspected origin of the lineations (Doyle, 1973, p. 5). 

3.2 The tracing of the Beech Grove lineament (Hollyday, Moore 
and Burchett, 1973) to the northern boundary of Madison County (Hollyday, 
personal communication) raises the question of its relationship to the 
trend of lineations (including the Anniston lineament of Drahovzal and 
others, 1974) in southwestern Madison County. A study of the bearings 
shows that they are at an angle to each other- -Beech Grove: N 20-25° E; 
Anniston: N 40-46° W. The Beech Grove trend, furthermore, enters 
Madison County near the Flint River and topographic analysis suggests 
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that It occury c-loci-licre lit Rid t son County. Lt, noucver, is not vell- 
developed, and is not easily identified (if at all (Resent) in the 
vicinity of the Anniston group of iineations. 

3.3 Reported rcrulto on studies in which remote sensing techniques 
are applied to subsurface hydrologic problems are rare. Coker, Marshall 
and Thompron (1969) working in limestone terrene in Florida, used thermo- 
graphy and color infrared photography and found that cooler subsurface 
temperature zones delineated areas of active sinkhole formation. These 
results suggest the ossibility that a similar effect might be expected 

as a result of otho subsurface processes involving ground water, and 
would act as a clue to subsurface conditions. 

3.4 The use of thermography as an indication of soil moisture 
has been addressed by many workers. For example, a recent report 
(Blanchard, Greely and Goettelman, 1974) contrasts the use of spectral 
reflectance and soil temperature methods to determine sell moisture and 
concludes that "the reflectance method appears promising for estimating 
soil moisture at the surface, and the temperature method appears prom- 
ising for estimating soil moisture near the surface (0 to 10 cm.) and 
perhaps slightly deeper." Either technique is applicable to the present 
research. 
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4.U Experiments CoiuJuc t •_■ < 1 During Phase '.two Research 

4.1 KH'-P Photographs 

An examination of photographs taken as part of the Skv lab 
Earth Resources Experimental Package (ERCP) was made. None proved to 
have Included the area of study. 

4.2 Hybrid Optical /Digital Data Processor 

An experiment to further minimize subjectivity in interpre- 
tation of linear features on ERTS or other imagery was conducted. The 
experiment involved the use of a hybrid optical/digital data processor 
in which a laser beam and horizontal bar filter were intended to produce 
a Fourier transform analysis. This is further explained in Appendix B - 
Lineation Detectimi System. 

Density contrasts in the imagery were intended to be trans- 
formed in such a manner that natural lineations were enhanced and 
printed by a computer-plotter combination, Scan lines on the imagery, 
however, were recognized by the optical data processor instead of natu- 
ral lineations, and the computer memory bank was overloaded. An optical 
method of de-emphasizing the scan lines may permit conversion of the 
ERTS imagery into a form compatible with the hybrid optical/digital 
data processor. 

4.3 Thermal Infrared Imagery 

4.3.1 Instrumentation 

The main thrust of the research was an experiment 
using thermal infrared imagery. Thermal radiation in the 8 to 14 micron 
range was recorded on 70 mm film by the Bendix LN 2 thermal scanner. 
Relative temperature differences between adjacent areas were evidenced 
by light and dark gray hues, cooler areas being darker, and warmer areas 
lighter. Soil moisture differences cause temperature differences because 
of water's capacity to retain heat. 

4.3.2 Rationale of the Experiment 

The rationale in using thermal Imagery to locate points 
of entry of potential contaminants of limestone aquifers is as follows: 



Cnt t*y of cont amt Hauls tlic l j* eyt«nie ,.>|<iilci‘ to more 

likely to oectti* where solution hoc occurroJ o« the liweatone bedrock 
surface and in the underlying, limeutonw. iStu-lt solution iwt often 
occurs along budding planes, or natural fractured in the rocky such as 
faults or joints 4 However i in nor h Alabama a rugolith (overburden , or 
"soil" in the engineering sense ) usually several lU’s of feet thick is 
derived from the insoluble residue loft as limestone or other carbonate 
rocks weathers. This overlies the bedrock and as a consequence solu' ion 
openings are hidden. 

To overcome the effect of the covering of the bedrock surface 
an indirect approach must be used. Textural differences in the soil pro- 
file which develops on the regol th probably occur and may be pronounced 
enough to permit such an indirect approach to the subsurface. The tex- 
tural differences develop in the following manner. In the weathering 
process, slightly acid water descends along fracture surfaces and removes 
the calcium carbonate comprising the limestone. This enlarges the frac- 
ture, allowing it to act as a drain. Then water tends to move 
toward the opening. The movement through the sands, silts and clays of 
the regolith causes graduul removal of some fine particles with a conse- 
quent increase in pore space. The greater pore space permits a greater 
volume of water to be present during a period of recharge after a rain. 
Conversely, it allows a lesser volume of water to be present after a 
period of uninterrupted drainage. In either case a moisture differen- 
tial would exist in the soil profile above a solution opening, as opposed 
to the profile above ait area where no solution opening is present. Inas- 
much as water has a heat holding capacity which can affect the soil temper- 
‘ature at the land surface, a moisture content differential may be reflected 
by a soil temperature differential that can be recorded by thermal infrared 
imagery. If thermal imagery is recorded under circumstances permitting 
the differences in temperature due to soil moisture to be shown, it becomes 
an aid in evaluation of subsurface geology and hydrology. Ground truth 
data is used to corroborate the information of the thermography. 



4.3.3 Availability of Tlu.i'nU) 1 'k .uuet* 

Hie themil scanner hecam* operatic in the late otvwee of 
the research period, eoitsequetU'ly , tlieriwgffipl^ was available for only 
otic season or meteorological condition - at the height of the growing 
season on June 20 and 27. Only the eastern half of the coverage was 
usable . 

4.3.4 Interpretation of Thermal Imagery 

Throughout most of the area, thermal differences recorded in 
in June 1974 arc attributable to the difference between the temperature 
of cropland and woodland. However, soil moisture differences may be 
responsible for short, linear thermal anomalies observable in areas 
judged to be covered with essentially homogeneous vegetation (Figure 1). 
The main value of the th 4 ;al features is in identifying the presence 
of presumed solution openings. This permits the taking of proper pre- 
cautions to prevent contact of potential contaminants with the land 
surface near the linear features . 

4.3.5 Ground Truth Data 

Ground truth in the form of soil temperature and moisture 
was recorded to help callibrate the imagery. Appendix A describes 
instrumentation and methods used in collecting the ground truth data. 

Figure 2 shows soil moisture determined gravimetrlcally by 
drying and weighing. Samples were taken a few hours after the over- 
flight of June 27. Molstyre ranged from 4.07. to 14.47,. Such moisture 
differences do not appear to be adequate to cause sufficient thermal 
differentials to override the thermal differences between crop and 
woodland . 

Figure 3 shows soil temperature determined by electrical 
resistivities recorded as the soil samples were collected for moisture 
determinations . 

4.3.6 Conclusions 

Under certain conditions (where vegetation is essentially 
homogeneous) temperature contrasts can be interpreted as indicating 
presumed fractures. By assuming greater solution of the limestone 



7 


bedrock ill t lie f frtCt»«VeO * yC«-<c al'fU'i coil be idunti I led where CuuL jii» 
iuuuts could find I iiuici* cut fy to the aquifer, ( bus pefKttt 1m* proper 
pi'ccauticHW to bu eXctcieeu. 

the i den t i flection of prcsw^O solutlonally-cnlorged ft*«cturce 
can also be used ay uti old in selecting locations foe U ft 1 liny to tap the 
ground water reservoir. 

5.0 Recommendations 

5.1 Experiments should be conducted in order to determine the 
optimum conditions under which to record the thermal imagery for the 
given purpose under investigation. Inasmuch as maximum thermal con- 
trast due to soil mo ire is desired, the imagery should be recorded 
under one or all of the following conditions: 

° Iti late autumn or early winter when vegetation is 
nearly dormant, so that the masking effect caused 
by the contrast between woodland and cropland is 
minimized , 

0 After a dry cold front lowering temperatures, for 
maximum temperature contrast between moist „nd 
dryer areas. 

° Two to three days after a widespread rain, for 
optimum difference in soil moisture. 

5.2 Red is sensitive to soil moisture occurring in red soils. 

Tills color soil occurs in north Alabama. (Consequently, photographs 
using natural color Ektachrome, taken when vegetation is dormant and 
•soil exposed to the maximum degree, should also be part of the experi- 
ment.) When exposed and processed so as to emphasize varying shades of 
red, the differences in soil moisture may be apparent. 

5.3 In order to test the supposition that differences in soil 
texture occur above solutionally-enlarged fractures in limestone, 
field studies should be made. Such studies should include locating, 

If possible, exposures of the bedrock surface which include solution- 
ally-enlarged joints. During field examination soil satnpj.es should be 
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collected vertically «K‘vu the solution leutuiVy u»vl 
them, It is i* ‘poctViJ I- 1 j.* • (• sot! it I*.\s v*J 1 I be in 
textural dif invitees not recor.“tw*«bJe In the Held, 
can be accomplished by soil auj.oriuy; on the r.roond o 
itieu occur In which aolotiounlly-eulacvcd t lectures 
from exposures to beneath the toil cover. 

Soil temperature and soil moisture data 

as needed . 


t a t ere 1 ly l re- < 
Jicatioe of any 
t’ini 1 ar sampling 
or face if opport ~ *■ 
can be projected 

should be collected 
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Append ic A 
CtVutTl Tcuth Pat a 
by U. E. McCartney 

Soil *■. ",>1cy were collected a few houm after the second overflight 
on June J7. The Jiten of collection were based on e .pec table locations 
of thermal (and moisture) anomalies an predicted from linoations on CRTS 
Imagery. 

Gravimetric methods cay used in determining soil moisture* The 
calculations Used weights before and after drying in the following 
formula : 

wet weight - dry weight 

Percent Moisture Content - ~ x 100 

ry weight 

Electrical resistivity cells were implanted near the collection 
sites a few days before collection. Measurements were made at the time 
of collection. Conversion of resistivity to temperature was possible 
through the use of callibration graphs. Several sequential determina- 
tions of soil moisture at a given site would allow preparation of 
correlation graphs of moisture versus resistivity (each correlation must 
be for a given site, because of differences in soil characteristics). 


Table A summarizes the data. 
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Table A 

Soil Moisture and I -•mpernture - Field Stations 







Soil 




Location 

Date 


Depth 

Temper- 

% of moisture 

Resistance 

Station 

(Sec-T-R) 

1974 

Hour 

(cm.) 

ature 

(by weight) 

(fl& 60°F) 

B 

SE-SW-SW- 

8-3S-2W 

6-27 

1100 

0-4 

n/a 

6.1% 

n/a 

C 

SW-SW-SE- 

13-3S-3W 

6-27 

1130 

0-4 

20°C 

8.47. 

l0 6.14 

D 

SW-NU-NW- 

31-3S-2W 

6-27 

1145 

0-4 

26 . 1°C 

9.87, 

1 0 6.72 

E 

NW-SE-Ntf- 

30-4S-2W 

6-27 

1200 

0-4 

24.4°C 

14.47. 

10 4.56 

F 

NW-NW-NE- 

5-5S-2W 

6-27 

1215 

0-4 

29 . 7°C 

4.57. 

10 6 

G 

SE-SE-NE- 

8-5S-2W 

6-27 

1240 

0-4 

25 . 6°C 

11.47. 

1q4 . 64 

H 

SE-NW-SW- 

16-4S-2W 

6-27 

1320 

0-4 

30 . 3°C 

10 . 5% 

1 & 6.17 

I 

SE-SE-NE- 
21 -4S-2M 

6-27 

1410 

0-4 

29.4°C 

10.07, 

1q6.44 

J 

NW-NW-NW- 

22-4S-2W 

6*27 

1425 

0-4 

n/a 

9.97, 

n/a 

K 

Field in which soil 

cell implanted plowed since 

installation 


L 

SW-NW-NE- 

14-4S-2W 

6-27 

1445 

0-4 

22 . 5°C 

9.27. 

10 6.35 

M 

NE-NE-NW 

35-4S-2W 

6-27 

1517 

0-4 

26.4°C 

9.07. 

io 5 ' 5 

N 

SW-NW-SW 

35-4S-2W 

6-27 

1530 

0-4 

25 . 8°C 

6.97. 

1q6.11 

0 

NW- NE-NE- 
3-5S-2W 

6-27 

1545 

0-4 

29.4°C 

4.07. 

io 6 - 24 

P 

SE-NE-NE- 

21-5S-2W 

6-27 

1600 

0-4 

30.3°C 

10.37, 

10 4.36 

Q 

SW-SW-NE- 

9-4S-2W 

6-27 

1355 

0-4 

n/a 

12.17, 

n/a 

R 

SW-Stf-NE- 
34-4S-2W ‘ 

6-27 

1613 

0-4 

n/a 

11.97, 

n/a 

S 

SE-SE-NE- 

36-4S-3W 

6-27 

1300 

0-4 

n/a 

13.37, 

n/a 


NOTE: n/a = not available (soil cell lost or not used) 
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Append ( >, u 

Lineal tun IVte-.-iion Svi . u ;, i 
by Joseph Kerr 

The* hybrid optical /digital prOOe.euor Used to attempt to identify, 
locate, and record lineal ions fi'ot>i CRTS imagery uus configured ay shown 
in Figure 4, 

The optical portion of the system performed the analysis necessary 
to identify 1 incat ions in the film. The control computer is a 11 P-ii 1 ltd) 
with a disc operating system supported by such peripherals as a Calcomp 
plotter, X, Y, Z modulated scope, A/D and D/A converters, and a digital 
voltmeter. Tills control computer scans the laser light across the sur- 
face of the film b controlling a set of X, Y galvanometer stages. The 
laser beam (1,5 mlllimoLer) is always perpendicular to the film plane 
because of the scanning laser optics whtcli is composed of the Galvoe, 
laser, parabolic mirror and the associated optics which are mounted on 
a fixture that rests on the Stable Table. 

The 1.5 millimeter beam of coherent light is diffracted by the ERTS 
image transparency which is located in the input film plane of the opti- 
cal processor. This diffracted light is gathered by the transform lens 
and is focused onto a rotational slit where its intensity is proportional 
to the magnitude of the Fourier transform of the image illuminated by the 
1.5 mm laser beam. The rotational slit is driven through 360° in discreet 
steps by a motor that is controlled by the control computer. Since linea- 
tions in the image appear as lines in the Fourier transform plane, the 
slit serves as a filter to locate the direction of any lineatioti that 
appears in the input film transparency. The photodiode provides an ana- 
log signal that indicates the magnitude of the light level. This analog 
signal is sent to the control computer storage disc via a digital volt- 
meter. The storage disc contains the position location, direction, and 
magnitude of each lineation detected. 



Al ter the entire i iltn t Cany pa reuey hau been scanned and analysed 
the cowtettty of the lineotion data on the disc arc a i sp 1 arod on the M 
Y, Z acopc so that the data pointy displayed are scaled to provide a 
one to one correspondence to the input irua&e . Those data cun also be 
plotted on the Cal comp plotter on a mylar- type material to provide an 
overlay that indicates the location of the 1 itieatiou with respect to 
the Image on the film transparency. 
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FIGURE 3 

SOIL TEMPERATURE (°C)ON27 JUNE 1974 
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FIGURE 2 

PERCENT SOIL MOISTURE ON 27 JUNE 1974 
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